o 

(N 



Generalization and applicability of the Landauer formula 
for non-equilibrium current in the presence of interactions 

H. Ness, x ' 2 '0L. K. Dash, 1 - 2 and R. W. Godby 1 ' 2 

department of Physics, University of York, Heslington, York YO10 5DD, UK 

2 European Theoretical Spectroscopy Facility (ETSF) 

(Dated: July 8, 2010) 

Using non-equilibrium Green's functions (NEGF), we calculate the current through an interact- 
ing region connected to non-interacting leads. The problem is reformulated in such a way that a 
Landauer-like term appears in the current as well as extra terms corresponding to non-equilibrium 
many-body effects. The interaction in the central region renormalizes not only the Green's func- 
tions but also the coupling at the contacts between the central region and the leads, allowing the 
total current to be further expressed as a generalized Landauer-like current formula. The general 
expression for the dynamical functional that renormalizes the contacts is provided. We analyze 
in detail under what circumstances Landauer-like approaches to the current, i.e. without contact 
renormalization, are valid for interacting electron-electron and/or electron-phonon systems. Numer- 
ical NEGF calculations are then performed for a model electron-phonon coupled system in order 
to validate our analytical approach. We show that the conductance for the off-resonant transport 
regime is adequately described by Landauer-like approach in the small- bias limit, while for the res- 
onant regime the Landauer-like approach results depart from the exact results even at small finite 
bias. The validity of applying a Landauer-like approach to inelastic electron tunneling spectroscopy 
is also studied in detail. 
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I. INTRODUCTION 

Electronic transport through nanoscale systems ex- 
hibits many important new features in comparison with 
conduction through macroscopic systems. In particular 
local interactions, such as Coulomb interactions between 
the electrons and scattering from localized atomic vibra- 
tions, become critically important. In crude terms, these 
effects are more important in nanoscale systems as the 
electronic probability density is concentrated in a small 
region of space; normal screening mechanisms are thus 
ineffective. 

It is most useful to have a simple expression for the 
electronic current or for the conductance of a nanoscale 
object connected to terminals. This is provided in the 
form an appealing intuitive physical picture by the Lan- 
dauer formula 1 , which describes the current in terms 
of the transmission coefficients of the central scatter- 
ing region and of distribution functions of the electrons 
in the terminals. However, in its original form the 
Landauer formula deals only with non-interacting elec- 
trons. This formalism has been used in conjunction with 
density-functional theory (DFT) calculations for realis- 
tic nanoscale systems 2 -"^ and has helped tremendously 
for the qualitative understanding of the transport prop- 
erties of such realistic systems. The apparent success 
of such approaches relies on the fact that DFT maps 
the many-electron interacting system onto an effective 
non-interacting single-particle Kohn-Sham Hamiltonian 
suited for the Landauer formalism for transport. How- 
ever when such a mapping becomes questionable for 
strongly-interacting electron systems, the original Lan- 
dauer approach has been found to be incomplete and 



unabl e to properly take into account the many-body 
effectPEl. 

The Landauer formula has been built upon by Meir 
and WingreerfSl to extend this formalism to a central 
scattering region containing interacting between parti- 
cles. It is then expressed in terms of non-equilibrium 
Green's functions and self-energies and in the most gen- 
eral cases it does not bear any formal rese mblan ce with 
the original Landauer formula for the currenP^ES. Other 
generalisation of Landauer-like approaches to include in- 
teractions and inelastic scattering have been developed, 
see for example Refs. |14ll6tH5] . 

It is therefore important to know the domain of validity 
of Landauer-like approaches in comparison to exact cur- 
rent calculations based on non-equilibrium Green's func- 
tions for treating electron transport through an interact- 
ing region connected to leads at different thermodynam- 
ical equilibria. This is the question we address in this 
paper by following a two step approach. 

First we reformulate Meir and Wingreen's work to once 
more express the current as the sum of a Landauer-like 
expression involving a transmission coefficient, plus a 
non-Landauer-like term arising from the non-equilibrium 
many-body effects. We further develop our theoretical 
framework to show that the interaction between parti- 
cles in the central region not only renormalizes the non- 
equilibrium Green's functions but also the coupling at the 
contacts between the central region and the leads. We 
hence obtain a generalized Landauer-like formula for the 
current in the same spirit as in Refs. |16I17I19I20J . How- 
ever our result for the dynamical functional that renor- 
malizes the coupling at the contacts is more general than 
the ansatz used in previous studies (Refs. 



Our result does not imply any constraints on the statis- 
tics of the non-equilibrium interacting central region. 

Secondly, we apply our theoretical framework to study 
a model system in the presence of electron-phonon 
(electron-vibron) interactions, connected to two non- 
interacting electron leads at non-equilibrium. We ana- 
lyze in detail the validity of Landauer-like approaches 
to describe the conductance and the inelastic electron 
tunneling spectroscopy (IETS) of such a non-equilibrium 
many-body interacting system. 

The paper is organized into two main sections. In the 
first (section [nj we develop our formalism to derive our 
generalized expression for the current. The implications 
of this are discussed in section |IIC| We then apply this 
formalism to the model system in section |IH| and show 
the results of numerical calculations. 



II. FORMALISM 
A. The model system 

Following Meir and WingreerJS], we consider a scat- 
tering central region (a quantum dot, a molecule, or a 
nanowire including interaction between particles) which 
is connected to two (left L and right R) leads. These 
leads are described by two non-interacting Fermi seas at 
their own equilibrium, characterized by two Fermi distri- 
butions /l(w) and fn(w). 

The Hamiltonian of the system is given by 



H= J2 £ a clc a + H int ({di},{d n };{al};{a x }) 



a=L,R 



^2 (Va,ncid n +K.C.), 



(1) 



n,a — L,R, 



where the summation indices a run over the left and right 
leads (L, R respectively) and depending on the choice 
of representation over momentum k or lattice site i in- 
dex, with $ a {c a ) creating (annihilating) anon-interacting 
electron, and {d^}; {d n } represent a complete, orthonor- 
mal set of states for the interacting electrons in the cen- 
tral region, and {a\}; {a\} represent a set of bosonic de- 
grees of freedom to which the electrons are coupled in 
the central region. These can be more or less extended 
phonons in a quantum dot or nanowire, or molecular vi- 
brations (vibrons) in molecules. 

There are two main approximations in the Meir and 
Wingreen approach to transport. The first is to consider 
that the interactions are localized within the central re- 
gion. This leads to specific properties for the self-energies 
used to calculate the non-equilibrium many-body Green's 
functions within the basis states of the central region 
only. The self-energies are then obtained as the sum 
of three contributions: two similar contributions arising 
from the electronic coupling of the central region to the 
left and right leads and the third arising from the inter- 
action between particles in the central region. 



The second approximation is to consider that the ini- 
tial correlations dies out in the long-time limit, and hence 
a steady state regime can be reached. It should be 
noted that a generalization going beyond the steady state 
regime has recently been given in Rcf. 22 . 



B. Non-equilibrium Green's functions and 
Landauer-like formula for the current 



In the steady state, the current II flowing at the 
left contact between the left lead and the central region 
is expressed in terms of three non-equilibrium Green's 
functions (the retarded G r , advanced G a and lesser 
G < Green's functions) of the dressed interacting central 
region^. 

Using the identity G> -G< = G r - G a and the defini- 
tion of the leads' self energies S^(w) = i/i(o;)ri(a;) and 
S^(w) = — i(l — /l(w))Fl(w), for which we recall that 
/l(w) is the Fermi distribution of the non-interacting 
left lead and Tl is obtained from the imaginary part of 
the retarded (advanced) self-energy E^ arising from 
the coupling of the central region to the left lead, i.e. 
r L (w) = T23mS2 /o (w), the current I L is given by 23 ^ 



lL = T J £ Tr[/LH FlH [Gr{uj) G<>)] (2) 
+ T L (u) G<(w)], 

where the trace runs over indexes n, m appropriately cho- 
sen to represent the electronic states of the central region. 
A similar expression can be obtained for the current Ir 
flowing at the right contact between the right lead and 
the central region by exchanging the subscript L -R- R 
in Eq. pi). For a current-conserving system, one then 
has I i, = —Ir. The famous result of Meir and Wingreen, 
(Eq. (6) in Ref. 13), is then obtained by evaluating the 
symmetrized current, / = (II — Ir)/2, to give 



/ = j I dw Tr [(f L (u)T L - f R {u)r R ) (G r (uj) - G») 



+ (r L H-r 7? M)G<H]. 



(3) 



Now, we can define more explicitly the specific property 
of the self-energies, namely additivity: E x (w) = E£(u;) + 
T, r (uj) + Tif nt (uj) where x is any component x — r,a, >, <, 
and the self-energies are defined within the central region 
by S^ (ui) for the coupling of the central region to the lead 
a, and by Ef nt (w) for the interaction between electrons or 
between electrons and phonons/vibrons. As mentioned 
in the previous section, the many-body interaction self- 
energy can be added to the leads' self-energies only be- 
cause the interactions are localized in the central region. 
Throughout the paper, we will also use a more compact 
notation for the leads' self-energy , i.e. E^+^j = E^ + E^.. 
Using the additivity property of the self-energy, and 
the fact that G<(w) = G r (w)E<(u;)G a (w) in the steady 



state regime, the symmetrized current / can be re- 
expressed as follows: 

I=j J~ {(fL-f R )Tr{T L G r r R G a ] 



Tr 
Tr 



j(y> _ y< \ 
(hT L ~ fRT R )G r l in * ind G a 



(T L - T R ) G 



r int r^a 



) 



(4) 



Introducing the non-equilibrium distribution function 
matrix for the interaction f^f defined from the inter- 
action self-energies as £< t = -f^f (£f nt - Sf nt ) (see 
Appendix [A]) , one can rewrite the symmetrized current 
/ as follows 

1 = T /£ ( (/l " /i?)Tr F^TrG*] 

T- \(fL-flf)T L -(f R -flf)T R ) (5) 



sir H^int ^int) pa 



G^iw) = [ g^^r 1 - E£ K ( W ) - S^M]- 1 . Depend- 
ing on the way the interactions are treated, the renormal- 
ization of the Green's functions may even go beyond the 
quasi-particle description of the interacting system. In 
any case, the important point is that T c ff already contains 
part of the electron-electron and/or electron-phonon in- 
elastic scattering processes. 

To complete our theoretical framework, we can make 
a further formal manipulation of the equations for the 
exact current, as given for example by Eq. pi, and end 
up with a more compact expression for the current which 
expresses a clear physical result: the interaction renor- 
malizes not only the Green's functions (G r,a ) but also the 
coupling at the contacts. 

To show this, it is more convenient to consider for the 
moment the current at only one contact (Jf, for example) , 
though one should not forget that in the steady state 
the current conservation implies II = —I R = I. The 
compact expression we find for II is the following: 



h = 



2e f dtu 

2tt 



(/lH - f R (u))Tr [T L G r T R G a ] (8) 



with the coupling to the right contact T R being renor- 
malized as 



The first term in Eq. p} and pi looks like a Landauer- 
like (LL) expression for the current, 



2e 



I hL = fJde(f L (e)-f R (e))T eS (e) 

= ^J^(fL-fn)Tr{T L G r r R G a ] 



(6) 



with an effective transmission 

T cff (e) = Tr [T L G r T R G a ] (e) = Trft (e)t(e)], (7) 

which can be interpreted with the intuitive physical pic- 
ture, as in the original Landauer formulation of electronic 
transport, in terms of transmission coefficients t(e) and 
propagation eigenchannels as defined in Refs. |37l38j . 

The second term in Eq. pi corresponds to non- 
equilibrium corrections due to the interaction. It is ex- 
pressed in terms of S in J. and of the different distribution 
functions. This term, not automatically small, cannot be 
recast in the form of extra transmission coefficients as in 
Refs. [37I38J . and already indicates in a way the break- 
down of the original Landauer formula for the current in 
the presence of interactional. 

One should note also that even if J LL looks like a Lan- 
dauer formula for the current with an effective transmis- 
sion T e ff(e), the interaction between particles is already 
taken into account in an exact calculation of the Green's 
functions. In this sense, 7 LL is not a conventional Lan- 
dauer current formula for single-particle elastic scatter- 
ing. The renormalization of the non-interacting refer- 
ence system is included in the retarded and advanced 
Green's functions via the corresponding self-energies: 



T fl (w) = r B (w)A(w), 


(9) 


W«)-C(w).„> 

f L {u) - f R {u) { lnt 


-£< t )H, (10) 



and 



A(w) 



where we recall that J(^(uj) is the non-equilibrium sta- 
tistical distribution for the many-body interactions as de- 
fined in Appendix [A) 

Equations pi and pfl0| (see also Eq. p3| in Ap- 
pendix [Bl represent the principal formal results of this 
paper. They imply that for an interacting central re- 
gion, one can always express the current in an general- 
ized Landauer-likc formula in which not only the retarded 
and advanced Green's functions are renormalized by the 
interaction but also the coupling at the contacts, as sim- 
ilarly found in Refs. [16I17I19T - I2T] . This generalized for- 
mula needs to be contrasted with the more conventional 
Landauer-like formula Eq. pi) in which the contacts of 
the central region with the leads are not renormalized by 
the interaction. 

Our expressions Eq. (8p0) are valid for any kind of 
interaction localized in the central region and general- 
ize the results of the previous studies (Refs. |16ll7ll9l - 
I21j ) because they do not imply any restrictions to the 
non-equilibrium statistics of the many-body interacting 
central region as we explain in detail in Appendix [Bj 

Finally, one recovers the more conventional Landauer- 
like formula (with no correction factors or equivalcntly 
with no renormalization of the contact couplings) when 
the quantity (£^ t (u;) — £^ t (u;)) vanishes, as can be 



clearly seen from Eqs. pi and Eqs. (8fl0). In the next 
section we discuss in detail the conditions for which 



this can happen. Note that the condition (E? t (w) — 
T,^ t (uj)) = does not necessarily imply that G r ' a (uj) = 
as well. Hence, even when the transport is well described 
by a Landauer-like formula ((E^ t (w) — E^ t (w)) = 0), nor- 
malization effects still occur and the transport is domi- 
nated by single-quasiparticle scattering. 



Discussion 



Clearly whenever (£^ t — £^ t )(w) = 0, there is no 
renormalization at the contact, and the current is sim- 
ply given by J LL . This may happen in two cases: ei- 
ther S ; ^ = for all w or only within finite range(s) of 



> 

int 



Ejj lt ) = should be included within the bias win- 



dow defined by the two Fermi levels \il and hr at non 
equilibrium. 

In order to understand how and why the quantity 



contour. 

In effect, any method which maps an interacting elec- 
tron/phonon system onto an effective one-particle (quasi- 
particle) scheme as for example in density-functional 
based technique (DFT, TDDFT) or other mean-field ap- 
proaches, will end up with no lesser and greater compo- 
nents for the corresponding self-energy describing the in- 
teraction. Hence a Landaucr-likc approach to the trans- 
port is entirely appropriate for such methods^'. However 
the mapping onto a one-particle scheme may not describe 
well strongly correlated electronic systems, A xc [n] being 
amenable to approximation, or is simply not possible in 
the general case of electron-phonon interaction. 

Now, let us turn to the second case: the interaction 
In the latter case, the relevant range of w for which spectral density SsmZ[ n t = j(£>. - £< t )/2 vanishes for 

one or more (connected) ranges of uj values. If this gap 
in $jm£[* nt is enclosed within the bias window then again, 
the current I will be determined only by the Landauer- 
like term 7 LL . 

As will be shown in detail below from numerical cal- 
culations for a model system, such a gap in QfrnSj"^ may 
exist in special cases of electron-phonon interaction. The 
gap in QfmEjJrt is then usually located around the Fermi 
level at equilibrium, and around the Fermi levels for the 
non-equilibrium cases at low applied bias only. These 
cases correspond to the regime studied by Imry et alW^ 
who derived a Landauer-like inelastic transmission for 
interacting electron-phonon systems and argued that the 
Landauer picture is still valid in the presence of interac- 
tion as long as multi-particle processes can be neglected. 

For electron-electron interactions, the situation is 
somewhat different. Even if collective excitations like 
plasmons present some qualitative bosonic analogy to 
phonons (there is a peak in the self-energy around the 
plasmon energy — as for the e-ph self-energy — and no 
much interaction spectral density elsewhere), there is al- 
ways however a non-zero contribution to the self-energy 
coming from the continuum of electron-hole excitations. 

The main difference from the electron-phonon in- 
teraction is that the phonon frequency loq imposes a 
restricted energy scale on the interaction, while for 
electron-electron interaction all energy scales are avail- 
able, making the corresponding interaction self-energy 
not vanishing, except in a infinitely small energy win- 
dow around the Fermi energy at equilibrium. So in prin- 
ciple, Landauer-like approaches for interacting electron 
systems are not valid for interactions treated beyond the 
mean- field/density- functional-based approximations^. 

Finally, one should note that the second term in Eq. ^ 
for the current, or the second term in the renormal- 



E ; „ t can vanish for an interacting system, let 
us first come back to the definition of the lesser and 
greater self-energies. These are specific components (pro- 
jections onto the real time axis) of the more general self- 
energy S int (r, t') with times r, r' defined on the Keldysh 
time-loop contour35E2HD Within the Keldysh approach, 
the lesser < (greater >) components of Ej nt (r,r') imply 
that the times t/t' are located on the forward/backward 
(backward/forward respectively) time-ordered branch. 
T, < ' > (u>) is simply the Fourier transform of E^^ (£,£') in 
the limit of the steady-state regime where any quantity 
depends only on the time difference X(t, t') = X(t — t'). 

First let us examine the first case: why would a self- 
energy have no lesser or greater components? For the 
so-called irregular self-energies^, we have the condition 
E(t, r') = E(r) S(t — t'). The self-energies for the in- 
teraction are instantaneous (local) in time. Hence they 
cannot have lesser or greater components, since the times 
have to be on the same time-loop branch. This con- 
dition of locality in time corresponds to two classes of 
physical effets. First when the self-energies describe one- 
particle potentials due to electron-electron or electron- 
phonon interaction, in other words it corresponds to the 
Hartree-Fock approximation for electron-electron inter- 
action and to only the Hartrec-like approximation for 
electron-phonon interaction. And second when the self- 
energies correspond to the so-called initial correlations 
which contain all contributions singular in time (see for 
example Refs. [43144] ) . 

There is also another class of problems for which there 
are no lesser or greater components of the self-energy. It 
is when the exchange and correlation effects for interact- 
ing electron systems are represented by an effective po- 
tential V xc (r, t) = 8A xc [n]/5n(r,t) being obtained from 
an exchange-correlation action functional A^n] of the 
electron density n(r, t). To this potential will correspond 
an effective self-energy that is local in both space and, 
more importantly, in timers* hence with no lesser and 
greater components for a generalization onto the Keldysh 



ization function A(w) in Eq. (10 1 involves the quantity 



(E? t — E ; „ t ) which is, in a series expansion of the in- 
teraction, proportional to the powers of the coupling 
constant (s) characterizing the interaction. In the limit 
of weak interactions, these terms represent small cor- 
rections to the Landauer-like current expression, and 
also give small contributions in the renormalization of 
both the Green's functions and the coupling at the con- 



tacts. Hence in the limit of weak interactions, con- 
ventional Landauer-like approaches could be confidently 
used and corrected by using perturbation theory for the 
interactiorPHUl 

Now we turn to presenting numerical calculations for 
a model system including electron-phonon interactions 
to illustrate our previous analysis. We compare results 
obtained from the exact current expression Eq. ^ or ^ 
with the current derived from the Landauer-like formula 
Eq. (pi), and by using different levels of approximations 
for the Green's functions. 



III. APPLICATION FOR INTERACTING 
ELECTRON-PHONON MODEL SYSTEMS 



In this section, we study in detail the validity of 
Landauer-like approaches for an interacting model sys- 
tem connected to two non-interacting electron reser- 
voirs at non-equilibrium. We concentrate on a model 
of electron-phonon interaction with the simplest version 
of the Hamiltonian (Eq. (fl])) for the central part: a sin- 
gle electron level coupled to a single vibration mode — 
the single-site single-mode (SSSM) modeP ^ whic h has 
also been considered in previous studies^ " 34 * 36 * . We 
also briefly describe below how to calculate the non- 
equilibrium Green's functions (NEGF) from this model 
Hamiltonian; the full theoretical details can be found 
elsewhere^!. 

We then apply our NEGF technique to the calcula- 
tions of the transport properties of the junction around 
equilibrium and out of equilibrium, and thus for differ- 
ent transport regimes. We analyze in detail if and when 
Landauer-like approaches can provide a good description 
of the transport properties in comparison to an exact 
calculation. 



A. Model Hamiltonian for electron-phonon 
coupling 

The Hamiltonian for the central region for the SSSM 
model is then given by 

He = eocTd + hujQa)a + 7o(a^ + a)d'd, (II) 

where one electronic level Sq and one vibration mode of 
energy w are coupled together via the coupling constant 
7o- 

Furthermore we choose a simple model for the struc- 
ture of the leads, which provides analytical results for 
the corresponding surface Green's functions, but in prin- 
ciple there is no particular restriction to be applied to 
the model or dimensionality of the leads. So in the fol- 
lowing, the left L and right R leads are described by 
two non-interacting one-dimensional semi-infinite tight- 



binding chains: 

— oo 
H L = ^ £ Lc\ c i + PL (c-Ci_i + C.C. 



+ ■00 



(12) 



h r = ^2 £ R c \ c i + Pr ( c ! c i-i + c - c -j ■ 



i=+l 



This model provides us with analytical expressions for 
the matrix elements of the leads' Green's functions at 
the terminal sites 50 : 



O) = e ik ^/P a , 



(13) 



with oj — e a + 2/3 a cosfc Q (u;), giving rise to semi-elliptic 
density of states of the terminal lead sites connected to 
the central region. 

The expression for the coupling of the central region 
to the L and R leads is then given by 



Vlc + V C r= Yj toa(4,d + rft 



(14) 



u=L,R 



with hopping integrals to a and c a= L = Ci—-i, c a= R 

Cl=+1- 



B. Non-equilibrium Green's functions for 
electron-phonon coupled system 

We use a non-equilibrium Green's function (NEGF) 
technique to calculate the propertie s of t he system in a 
similar manner to previous studiea^ " 34 * 36 l The details of 
our NEGF calculations are described in detail in Ref. [49] , 
but we briefly summarize our application of them here. 

The Green's functions are calculated via Dyson-like 
equations for the retarded and advanced Green's func- 
tions G r ' a (w): 

G^ a {u) = <£» + .g^HE^HG^H, (15) 

where g^ 11 ' is the non-interacting Green's function for the 
isolated central region. 

For the greater G > (lu) and lesser G < (u>) Green's func- 
tions, we use a quantum kinetic equation of the form 



G° 



(I + G r S r )ffc I 1 + SQGa ) + G r E > - < G a . (16) 



Here T, x (oj),(x — r, a, >,<) is a total self-energy con- 
sisting of a sum of the self-energies from the constituent 
parts of the system: 



E*( W ) = E|( W ) + ES(«) + 5^ t («). 



(17) 



The leads' self-energies Ef, ^.(w) = E|(w) + T, R (ui) aris- 
ing from the non-interacting leads a = L, R are given 
by 

Z<(w) = -2i$Sm[Tr a {u)]f a (u), (18) 

E>(«) = -2i9top;( w )](/ a (u,)_i), 



where gQ a is given by Eq. ( 
distribution for lead a, wit 



13 ) and f a is the Fermi-Dirac 
eq. 



:h Fermi level fx a = /j, cq + ri a eV 
and temperature T a . At equilibrium, the whole system 
has a single and well defined Fermi level fi cq . Out of equi- 
librium, a finite bias is applied throughout the junctions. 
Within our model Hamiltonian, the fraction of electro- 
static potential drop at the left contact is t/l = ±rjy 
and r/R = +(1 — rjy) at the right contact, hence 
Vl —VB^eV is indeed the applied bias, and r\y € [0, 1]. 
The self-energy for the interaction in the central re- 
gion, £ int (w), is obtained from a non-equilibrium many- 
body perturbation expansion^ of the electron-phonon 
coupling term in the Hamiltonian, Eq. (11). As for 



a conventional many-body perturbation expansion, the 
self-energy is associated with a series of Feynman di- 
agrams for the interaction. In the current work, we 
consider only the lowest order diagrams, i.e. the Born 
approximatio n (BA) or equivale ntly the Hartree-Fock 
approximatio n l 23 | 24 | : ^ 1 281 |34 | 36 | 5 2i53j rp Qe exac ^ expressions 

for Ef nt , (x — r, a, >, <) at the Hartree-Fock level and be- 
yond are given in Ref. [49] and we do not reproduce them 
here. 



C. Numerical results 

We divide the calculations into two different types 
of transport regimes. The first of these is when ei- 
ther £o <^ fi cq or £o 3> A i ° q - known as the off-resonant 
regime. It corresponds to a poorly conducting junction 
(i.e. semiconductor or insulator-like behavior) dominated 
by strong tunneling at low bias. The second transport 
regime is when Eq ~ /i cq and known as the resonant 
transport regime. This regime corresponds to a good, 
metallic-like, conducting junction. 

We will see below that depending on the nature of the 
transport regime, the Landauer-like approaches may be 
sufficient, under certain conditions, to describe the con- 
ductance G(V) = dl/dV or the inelastic electronic tun- 
neling spectroscopy (IETS) properties of the junctions. 



The self-energy quantity (E? t — EJ^ t ) 



As discussed in section II C the quantity (£• 



£< t ) 
plays the key role in determining whether or not the 

Landauer-like approaches are valid. This quantity is plot- 
ted in Figure [l] for the off-resonant transport regime and 
for intermediate electron-phonon coupling strength. Cal- 
culations were performed self-consistently using the low- 
est order electron-phonon diagrams, i.e. within the con- 
ventional self-consistent Born approximation (SCBA) 49 . 
By definition £> t - £< t = E[ nt - Sf nt is a purelyimag- 
inary function for conserving approximations^^!], it 
presents features (peaks) corresponding to the excita- 
tions of the system. The features obtained for non- 
equilibrium conditions, especially when real excitations 
can be created in the system (applied bias V > ujq), are 



A 




FIG. 1: The (E> t - E£ t )(w) quantity for the SSSM model 
within NEGF-SCBA for off-resonant transport regime and 
intermediate electron-phonon coupling strength 70 /wo = 0.65. 
The curves are obtained for different applied bias eV = \il — 
fj,R, and are offset vertically for clarity. The applied bias is 
also shown and given by the chemical potentials of the left 
jj,l (left-pointing arrows) and right (ir (right-pointing arrows) 
leads respectively. The other parameters are so = +0.5, ujq = 
0.4,7o = 0.26,£ol„r = 0.2, T L , R = 0.011, r) = 0.025, rjv = 1. 
At zero and low bias (here V < 0.25) the interaction self- 
energies difference (E? t — E^ t ) is zero within the bias window, 
and hence the current is given by a Landauer-like formula. 



strongly different than the features obtained for equilib- 
rium (no applied bias V = 0). 

At zero and low bias (V < 0.6 loq for the set of param- 
eters used in Figure [lj, the difference between the inter- 
action self-energies (E{^ t — E^ t ) is virtually zero within 
the bias window \pu, hl\- This means that the current in 
Eq. (|5| is effectively given only by the first term J LL , and 
hence Landauer-like approaches are sufficient to describe 
the transport properties of the system for such biases. 



2. Conductance and inelastic electron tunneling 
spectroscopy (IETS) 

In order to analyze in detail the different contribu- 
tions to the conductance and the conditions for which 
the Landauer-like approaches can be valid, we have per- 
formed calculations for the current within different lev- 
els of approximation. In the following, we consider four 
different kinds of approximation: firstly, in the absence 
of interaction, the current is obtained from the non- 
interacting Green's functions G r Q a {u) = [g^ a (uj)^ 1 — 
^2+ij( w )] _1 ana - corresponds to the original Landauer 
formulation 7 LL [Go] calculated with Eq. (pi). 

Secondly, the Landauer current I LL [Go] can be cor- 
rected to include the interaction effects at the lowest 
order of the coupling parameters, as in perturbation 
theory extended to non-equilibrium conditions. This is 
done by calculating the current in Eq. p} using only 
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FIG. 2: Dynamical conductance dl/dV from the exact expres- 
sion for the current I(V) and the corresponding Landauer-like 
current I LL (V). Green's functions calculations are performed 
self-consistently for the electron-vibron coupling within the 
Born approximation and for intermediate coupling 70 /wo = 
0.65. In the off-resonant case (upper panel, eo — +0.5), there 
is a good agreement between the two conductances for the lin- 
ear regime and for low biases, i.e. when (E^ lt — E^ t )(w) ~ 
as shown on figure fT] Hence Landauer-like approaches are 
valid in the low bias regime for off-resonant transport. For 
larger biases (V > 0.5), dI Lh /dV gives unphysical negative 
conductance values. In the quasi-resonant case (lower panel, 
£0 = +0.2), only the linear conductance is well reproduced 
by 7 LL . Large deviations of d/ LL jdV from the exact result 
occur at small applied bias. Hence the Landauer-like fornrula 
for non-equilibrium current is not valid for the quasi-resonant 
case. The other parameters for the calculations are the same 
as used in figure [T] Wo = 0.4, 70 = 0.26, toL,R = 0.2, -qv = 1. 
The insets show the ratio of the conductances calculated with 
the exact and the corresponding Landauer-like current. 



the non-interacting Green's functions Gq" in the first 
term of Eq. Q and in the evaluation of the interac- 
tion self-energies SJ^j. [Go]; we denote this current by 

^perturb = I [Go] + A/[£i n t [Go]] . 

Then the last two kinds of approximations include full 
renormalization effects in the Green's functions. The first 
of these corresponds to a non self-consistent BA calcu- 



lation using the self-energies E; nt [Go] to renormalize the 
Green's functions as follows: G r ^ a A {uj) = [Gq"^) -1 — 
E£j£[Go(w)]] -1 . These Green's functions are then used to 
calculate the current I[G B a] = / LL [G BA ] + A/[S lnt [G ]]. 
Finally, the last approximation, correspond to a fully self- 
consistent renormalization SCBA calculation performed 
as described in Ref. [15] , from which we obtain the exact 
current I[G S cba] = I L ' L [G S cba\ + AI[T, int [G SC BA]}- 

The dynamical conductance G{V) is obtained as usual 
from the first derivative of the current versus the applied 
bias G(V) — dl/dV. Typical examples for both the off- 
resonant and resonant transport regimes are shown in 
Figure [2] where we compare the conductance obtained 
from the exact current with the Landauer-like current 
I LL (V) using full SCBA calculations. 

For the off-resonant transport regime (Fig. [2f a)), 
which is dominated by strong tunneling at low biases, 
there is a good agreement between the exact conductance 
and the Landauer-like conductance for both the linear 
regime and the non-linear regime at low biases, i.e. when 
(S? t — E^ t )(w) ~ (as shown on Figure nj). Hence in 
this case, Landauer-like approaches are valid to describe 
the tunneling transport properties in the low-bias regime 
for off-resonant transport. However, strong discrepancies 
between the two conductances occur for biases around 
the first renormalized electronic resonance V ~ 0.3, even 
before real excitations of the phonon mode are available. 
For even larger biases, (V > 0.5), dI Lh /dV gives unphys- 
ical negative conductance values. 

For the quasi-resonant transport regime (Fig. ^[b)), 
where the renormalized electronic resonance eo— 7o/wo ~ 
is close to the Fermi level at equilibrium, only the linear 
conductance is well reproduced by I LL . This is essentially 
due to the fact that within conservation approximations 
for the sclf-cncrgics, the linear conductance is not renor- 
malized by the interaction 5 ^ 5 -^, as we have also shown in 
detail in the appendix of Ref. [49] . However large devia- 
tions of the Landauer conductance from the exact result 
occur quickly at small applied bias. Hence Landauer-like 
approaches for non-equilibrium current are not valid for 
the quasi-resonant case, and probably not as well for the 
resonant case. 

We now turn to the inelastic electron tunneling spec- 
trum (IETS) which gives information about the selective 
excitation of the system. The IETS is usually obtained 
from the second derivative of the current with respect 
to the applied bias d 2 I/dV 2 . The IETS curves present 
features, peaks or dips^S, at biases corresponding to the 
energy of a specific excitation, in our case to the energy 
of one or several excitations of the vibration mode mo. 
Being the derivative of the conductance, the IETS curves 
also present features at biases corresponding to peaks in 
the conductance (see for example Figure [2]). 

We have found that in order to get a better aspect ra- 
tio for the IETS features corresponding to phonon excita- 
tions, it is more convenient to normalize the IETS curves 
by the dynamical conductance, i.e. [d 2 1 / dV 2 } / [dl / dV] = 
d/dVlnG(V). Typical examples for both the off- 
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FIG. 3: IETS signal d 2 I/dV 2 , normalized by dynamical con- 
ductance G(V) = dl/dV, obtained from the exact expres- 
sion for the current I(V) and the corresponding Landauer-like 
term I hh (V). Calculations are performed for the off-resonant 
transport regime and the strong coupling limit 70/^0 = 0.875 
using different level of approximations for the Green's func- 
tions, as explained in the legend (see main text for detail). 
The inelastic vibron excitations are present in the IETS de- 
rived from the exact I(V) and are located at integer mul- 
tiples of the vibron energy wo- They give a positive con- 
tribution to the baseline, as expected for opaque junctions 
(insulator/semiconductor-like behavior). However this signal 
is absent in the IETS derived from the Landauer-like term 
7 LL (y). This implies that even with normalized Green's func- 
tions the Landauer-approaches are not able to correctly repro- 
duce the inelastic features in IETS for the off-resonant regime. 
The parameters for the calculations are £0 = +1.5, Wo = 
0.4, 70 = 0.35, t 0L ,R = 0.11, T LM = 0.011, jj = 0.02, rjv = 1. 



resonant and resonant transport regimes are shown in 
Figures [3] and [4] respectively. 

We compare the IETS signals obtained from the four 
approximations used to calculate the current: ^Landauer 
for the Landauer current of the non-interacting system, 
^perturb for the Landauer current corrected by first order 
perturbation theory for the electron-phonon coupling, 
and full rcnormalization within a non-self-consistent and 
self-consistent scheme I[Gba] and I[Gscba] from which 
the corresponding Landauer-like contribution can be ex- 
tracted. 

We first comment on the results obtained for the off- 
resonant case, Figure W[&)- As expected, the original 
Landauer approach does not provide any feature at the 
bias V = ujq since there is no interaction. The IETS 
signal obtained from the fully self-consistent calculations 
shows however two features (peak s in the case of strong 
tunneling regime at low bia d 25 * 52 *) at biases V = ujq and 
V = 2ojq. They correspond to inelastic processes in- 
volving the excitation of one and two vibration modes 
respectively. 

The rising background of the curves for bias V > 1.2 
corresponds to the feature associated with the main res- 



onance in the conductance (see the corresponding main 
central peak in the conductance curves in Figure pi). The 
IETS curve rises at lower bias for the exact calculation 
compared to the calculations for the non-interacting sys- 
tem, simply because the exact calculations include a full 
renormalization of the electronic level £o- Such a renor- 
malized level is then shifted towards lower energy by the 
full dynamical polaron shift 49 60 . 

It is interesting to note that the first-order perturba- 
tion correction to the Landauer current, given by the 
second term in Eq. J5 1 when evaluated from Go only, pro- 
vides not only a qua itatively good feature in the IETS 
signal at V = ujq but also a partial renormalization of the 
electron resonance. That is, the background of the curve 
rises faster than for the non-interacting case, and this 
corresponds to a shift of the electron resonance towards 
lower energy by a partial polaron shift. 

Calculations performed non self-consistently provide a 
partial renormalization of the electron resonance, how- 
ever this is closer to the exact result than that obtained 
from perturbation theory. Additionally, the correspond- 
ing IETS signal also shows only a feature at V — cj as 
for perturbation theory, but its general aspect is again 
closer to the exact result. 

Another important result of our calculations is that the 
IETS signal calculated from only the Landauer-like term 
in the current I[Gba] and I[G$cba] does not contain 
any features at V = loq or V = 2wo, as shown in Figure 
|3j Although the IETS has the correct rising background 
and shows the correct corresponding renormalization of 
the electronic level, it fails to correctly reproduce the 
features associated with inelastic processes. Hence in the 
off-resonant transport regime the renormalization (self- 
consistent or not) of the Green's functions G ra , from 
which the Landauer-like current is obtained, does not 
contain the appropriate physical information to correctly 
describe the corresponding IETS signal. However, as we 
have seen above, it is good enough to describe the overall 
behavior of the conductance at low biases. 

Now we turn to the resonant transport regime, and 
check if the trends obtained for the off-resonant regime 
hold here as well. It should be noted that in the following 
calculations, we have considered the resonant transport 
regime in the case of strong coupling of the central region 
to the leads, toL,R ~ Pl,r- We are then dealing with an 
almost homogeneous one-dimensional system with metal- 
lic like behavior at equilibrium, in which the propagat- 
ing electrons are coupled locally to a single localized vi- 
bration mode. We have chosen this somewhat peculiar 
regime so that the features in the IETS signal associ- 
ated with the inelastic processes are not 'distorted' by 
the features associated with resonant-like transport. In 
other words, the background of the IETS signal around 
the excitation energies is fairly fiat. 

The corresponding IETS curves are shown in Figure HI 
for the different kinds of approximation used to calculate 
the current. As expected, the IETS signal calculated for 
the non-interacting system does not show any feature at 
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FIG. 4: IETS signal d 2 I/dV 2 , normalized by G(V), obtained 
from the exact current I(V) and from the corresponding 
Landauer-like 7 LL (V) term. Calculations are for the reso- 
nant transport regime, intermediate electron-vibron coupling 
7o/cjo = 0.65 and strong coupling to the leads. The different 
approximations used to calculate the Green's functions are 
shown in the legend (see main text for detail). The inelastic 
vibron excitation is present in the IETS signal derived from 
the exact I(V) and is located around the vibron energy ojo. 
It corresponds to a negative contribution to the baseline, as 
expected for mostly transparent junctions (metallic-like be- 
havior). Interestingly, this feature is also present in the IETS 
derived from 7ll(V) in contrast to what is obtained for the 
off-resonant case. Hence, for resonant transport, it seems that 
Landauer-like approach can reproduce the inelastic IETS fea- 
tures at V = ojo . The parameters for the calculations are 
eo = 0,u;o = 0.3,7o = 0.195, t 0LyR = 1.50, T L , R = 0.011,7? = 
0.025, Tjv = 1. 



FIG. 5: IETS signal d 2 I/dV 2 , normalized by G(V), obtained 
from the exact current expression I(V) and from the corre- 
sponding Landauer-like I hh (V) term. The Green's function 
calculations are performed self-consistently for the resonant 
case and stronger coupling (70/cJo = 0.875) than in Figure El 
Unexpectedly, even for such a strong coupling, there is no fea- 
ture in the IETS signal at V = 2ojo obtained from the exact 
current. The Landauer-like current provides a small feature 
at V — 2ujo with however a spurious positive contribution to 
the baseline. The other parameters are eo = 0,ujo — 0.4, 70 = 
0.35, t 0L , R = 1.40, T LM = 0.015, 77 = 0.017, r]v = 1. 



scribe the IETS signature of the inelastic process at the 
lowest excitation energy. 



V = ujo, while the IETS signal obtained from fully self- 
consistent calculations show a dip at V = Uq. Such a neg- 
ative contribution to the base line is to be expected in the 
case of good conductor d 24 * 25 * for which electron-phonon 
coupling is associated with electron backscattering. The 
results obtained from non-self-consistent renormalization 
are very similar to the exact results. Interestingly, the re- 
sult obtained from perturbation theory gives a feature in 
the IETS signal at the right bias, but however with the 
wrong sign. 

It appears that in all the cases we have studied, first 
order perturbation theory always gives a positive con- 
tribution (i.e. a peak) to the IETS which is generally 
incorrect. Indeed it has been shown that th e inelastic 
features of the IETS can be both peaks or dip d 24 | 25 | 36 | 6 1J 
depending on the nature of the conductor and essentially 
on all the parameters describing the systerrPn 

The most interesting result shown in Figure [4] is that 
the IETS signal obtained from only the Landauer-like 
term in the current I[Gba] and I[Gscba] shows the ap- 
propriate feature at at V = ojq. Hence for the resonant 
transport regime with strong coupling to the leads, the 
renormalization of the Green's functions from which the 
Landauer-like current is derived is good enough to de- 



Finally, we want to comment on the features in the 
IETS observed at V = 2ojq, which should correspond to 
the excitation of two vibron modes by the injected non- 
equilibrium charge carriers. Such a feature is observed in 
the rising background of the IETS in the strong tunneling 
regime (off-resonant transport regime) (Fig.pl; however 
this feature is clearly absent in the resonant transport 
case. Even in the case of strong electron-vibron cou- 
pling, shown in Figure [5j there is basically no feature at 
V = 2ojq in the IETS signal derived from the exact ex- 
pression of the current, although there is a small feature 
in the IETS signal derived from the Landauer-like current 
expression. The absence of such a feature at V = 2w for 
resonant transport is in agreement with previous studies 
(Refs. [24162] ). However at the moment there is no sat- 
isfactory physical explanation for the absence of such a 
feature at V = 2cuo. We postulate that such an absence 
may be due to the partial resummation of the electron- 
phonon diagrams in the many-body Green's functions. 
As we have already shown in Ref. |3S], higher-order di- 
agrams for the interaction (beyond Hartree-Fock / Born 
approximation) play a very important role in correctly 
describing the properties of the electron-phonon coupled 
system, even at intermediate electron-phonon couplings. 
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IV. CONCLUSIONS 



Acknowledgments 



In this paper, we have re-addressed the problem of 
the breakdown of Landauer-like approaches for electronic 
transport in the presence of many-body interactions. 
Starting from the original work of Meir and Wingreen^ 3 -', 
we have once more expressed the current as the sum 
of a Landauer-like expression (I LL ) involving the con- 
cept of single-particle transmission probabilities plus a 
non-Landauer-like term arising from the non-equilibrium 
many-body effects. We have further developed our the- 
oretical framework to show that the interaction in the 
central scattering region renormalizes not only the non- 
equilibrium Green's functions but also the coupling at the 
contacts between the central region and the leads. We 
have hence obtained a new form of generalized Landauer- 
like formula for the current, Eqs.(8fl0), in a similar way 
to Refs. |16|17|19|20) . However our result for the dy- 
namical functional that renormalizes the coupling at the 
contacts is more general than the Ng ansatz2B used in 
Refs. |16ll7ll9l2Tj] . Moreover our result does not impose 
any constraints on the statistical properties of the non- 
equilibrium interacting central region. 

We have then applied our theoretical framework to a 
model system of electron-vibron interacting nanojunc- 
tion. We have analyzed in detail the domain of validity 
of Landauer-like approaches, i.e. without renormaliza- 
tion of the contacts, to describe the conductance and the 
inelastic electron tunneling spectroscopy (IETS) of such 
a non-equilibrium many-body interacting system. 

Our results confirm that generally Landauer-like ap- 
proaches are not adequate to describe the transport prop- 
erties of such interacting systems for the whole range 
of applied biases (linear to highly non-linear regime) 
and for all the transport regimes (good, metallic-like 
to mediocre, insulator-like conductors). In general, the 
correct transport properties are only obtained from ex- 
act non-equilibrium many-body Green's function calcu- 
lations. However, there exist a certain number of condi- 
tions in which a Landauer-like approach can reproduce 
fairly well either the conductance or the IETS signal. 
For example, renormalization of the Green's functions 
in a Landauer-like approach is already sufficient to ac- 
count qualitatively for the inelastic features in the IETS 
signal for the resonant transport regime, i.e. / LL gives 
features at ujq in the IETS signal. However this is not the 
case for the off-resonant transport regime, for which the 
Landauer-like current J LL fails to reproduce the inelastic 
features in the IETS signal. 

Finally, we believe that higher order diagrams, as stud- 
ied in Ref. [49], may change the detailed features of 
the IETS signal, especially for higher energy excitations. 
However the accuracy to/with which the Green's func- 
tions are calculated (Born approximation/Hartree-Fock 
or beyond) does not alter the main conclusions of our 
work concerning the applicability of Landauer-like cur- 
rent formula versus exact derivation of the current. 
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Appendix A: Relationships between Green's 
functions 



In order to keep the formalism simple, let us first con- 
sider that the Green's functions and the corresponding 
self-energies are simply complex functions of w, i.e. we 
are dealing with an interacting central region containing 
only one site/one electronic level (SSSM model). 

When the system is at equilibrium (/& = fn = / eq ), 
the lesser (greater) Green's function is related to the ad- 
vanced and retarded Green's functions : 

G<< oq H = -/ cq H {G r < cq (uj) - G a ' cq (oj)) , (Al) 

and 

G>' cq H = -(/ oq (w)-l) (G r ' cq (w) - G a < oq H) . (A2) 

These relationships are at the center of the fluctuation- 
dissipation theorem for equilibrium. They can also be 
recast as follows, 



G>< cq H = ~ e ^-^)/ kT G<< cq (u;), 



(A3) 



and they then define a relationship between the greater 
and lesser Green's functions for statistical averages 
in the grand canonical ensemble at finite tempera- 
ture (the so-c alled Kubo-Martin-Schwinger boundary 
condition d 35 * 54 * ) . Similar relationships also exist for the 
self-energies £<>>> r > a (see for example Ref. [63]). 

For non-equilibrium conditions, there is no unique 
Fermi level at finite bias (or no unique temperature if 
Tj, 7^ Tr) in the whole system, and the relationships 

A3[) no longer hold. This is 



given above by Eqs. (Al 



an important feature of the non-equilibrium formalism 
for which conventional equilibrium statistics need to be 
reformulated^. However, the self-consistent calculations 
of the Green's functions and self-energies in the non- 
equilibrium case permit us to define new non-equilibrium 
distributions. For example, the non-equilibrium distri- 
bution / NE (w) is defined from the Green's functions as 
follows: 



G<{uj) = -/ W »(G» - G»). 



(A4) 



This definition is reminiscent of the so-called Kadanoff- 
Baym ansatz which has been generalized to the non- 
equilibrium conditions and to the time representation of 
the Green's functions (see for example Refs. [65166] ). 
Similarly, we also define the non-equilibrium distribution 
f^f(uj) from the interaction self-energies as follows: 






] int(^) : 



-/£?(«)(££*(«) - 



=£*(«))■ (A5) 
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There is no a priori reason for these two non- 
equilibrium distribution functions to be equal to each 
other at non-equilibrium. Both distribution functions 
contain information about both the non-equilibriu m and 
the many-body interaction effects in the systerrP2E3 
However, at equilibrium, these distribution functions are 
as expected equal to each other and to the conventional 
Fermi-Dirac equilibrium statistics / NE = /^ t E = f cq . 

As an example, the non-equilibrium distribution func- 
tion /q E (uj) for a non- interacting system^ is given by 
the weighted averaged of /l,r by the coupling at each 
contact Tlr: 

,ne/, a _ /lMFlM + f R (u})T R (uj) 

Jo K 10 ) — ^ I \ . t-i i \ ' v °' 

T l (uj) +T r (uj) 

The asymptotic values of distribution functions are de- 
fined from behavior of the Green's functions and self- 
energies at large w, and follow the conventional statistics: 



pNE 



NE/ 



r>) = /&» = /o » = /«n = n(u) 




(A7) 



At equilibrium, one recovers the equilibrium statistics for 
any distribution function Jq E = / cq . 

Now, we can generalize our formalism to central re- 
gions containing several electronic stat es. The distribu- 
tion functions then become matrices f l 13 l 45 l 69 l with ele- 
ments f nm given by 



A„, 



i 



fnl \X ln 






(A8) 



where X is either a Green's function G or a self-energy 
E, and the indices n,m are appropriate indices to label 
the electronic states of the central region. 



s^r L 



After noticing that E^T^ 

iE^ E|~ , R , we can finally obtain a compact form for A(o 



iS^S>, 



A(w) 



E<HE>H 



s>he<h 



;ME£ +fl (w)-E>( W )E< ( w ) 



(B3) 



which is another way of expressing the important result 
of this paper given in Eq. ( 10 1. 



Now we are going to relate our principal results to 
previous studies using the Ng ansatJSE!] -j^g ]\jg 
ansatz, developed to study the Anderson model out of 
cquilibriunPSHED j s based on using an apparently more 
convenient way to express the full lesser (greater) self- 
energy in terms of the lesser (greater) self-energy for the 
non-interacting system: 



'(«) 



Stfr^w) A(w) 



where Ti x = £ 



L+R 






J L+R\ 



and E 



(B4) 



L+R 



— ^L 



J /V- 



and 



A is a dynamical "renormalization" quantity, to be de- 
termined from the condition E> — H < = E r — E a . We 
will show below that this ansatz actually implies strong 
constraints on the statistics on the non-equilibrium in- 
teracting systems. 

By using the Ng ansatz to expr ess the lesser and 
greater self-energies E <:> (w) in Eq.(B3l, one can eas- 
ily see t hat our renormalization functional A(w) given by 
Eq.(|B3j) is just equal to the dynamical quantity A(w) 
of the Ng ansatz: A(w) = A(w). We then recover 
all the expressions for the current previously derived in 
Refs. |16I17I19W2T] from our main results Eqs. (8fT0l. 



However, the Ng ansatz presents some intrinsic lim- 
itations. To prove this, it is sufficient to calculate the 
non-equilibrium distribution function /^ E within the Ng 

NE/ 



ansatz. Starting from the definition of /^ E (u;), i.e 
/* E M = -E< t ( w )/(Er nt H - Ef nt M) , 



(B5) 



Appendix B: The Ng ansatz 

Here again, we consider in the following mathematical 
developments that the Green's functions and the self- 
energies are simply complex functions. Extension to ma- 
trices is rather straightforward but must be done with 
care using the notations and definitions given in the main 
text and in Appendix |A") 

Using the definition of f^f(u>) ■ 



E< t H = -/&» (elm - E? nt M) 



/^H(E> t H-E< t H), 



(Bl) 



and the fact that f a {oj) = — iE^(uj)/Tl(ui), we find af- 
ter more formal manipulations that the renormalization 
functional A(w) (Eq. ( 10 )) can be re-expressed as: 



A(w) 



; ^M£> t 



M 



s£M££ t 



(«) 



E<( W )r iJ (o;) 



E<( w )r L ( w ) 



(B2) 



it is straightforward to show that f^ is then given by 



/£?(«) 



E< +Ji ( W )(A(u;) - 1) 



\^L+R. '-'LA 



-E< +i? )(A-l) 

/i(«)r z (w) + / fl (w)r R (w) 



(B6) 



T l {uj)+T r (oj) 



= /<TM 



the non-equilibrium distribution function for the non- 
interacting system! 

Just to confirm the consistency of our derivations, if 
we use the above result J^ e (uj) — Jq E (uj) in the defi- 
nition of our renormalization functional A(w) given by 



Eq. (10) and the Ng ansatz for E>f , i.e. E>^(w) = 



E^YTf^w) (A(w) — 1), we find again and consistently that 
A(w) = 1 + (A(w) — 1) = A(w), as expected. 

However we have founcP^that the condition f^ E (uj) — 
/ NE (w) implies necessarily that / NE (w) = / NE (u;). In 
other terms, the Ng ansatz implies that the full non- 
equilibrium distribution / NE of the interacting system, as 
well as /^[ E , are equal to the non-equilibrium distribution 
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function for the non-interacting system! This is a condi- 
tion that is in contradiction with the fact that both distri- 
bution functions should simultaneously include both the 
non-equilibrium effects and the many-body interaction 
effects. In fact f NE (w) is actually a functional of both 
the non-equilibrium distribution function for the non- 
interacting fQ E (uj) and the many-body interactio n 1 64 * 67 !. 
Hence we conclude that our expression for the renor- 
malization of the coupling at the contact A(w) given 
by Eq. ( 10 1 is more general than the definition used in 
the Ng ansatz. The latter is not taking fully into ac- 
count the interaction effects at non equilibrium. It actu- 
ally corresponds to a lowest-order expansion of the full 
non-equilibrium distribution in terms of only the non- 
equilibrium distribution function of the non-interacting 
systerrp2. Or in other words, the Ng ansatz consid- 
ers that the statistics of the interacting central region 
is dominated by that of the non-interacting leads at non- 



equilibrium, and that the interaction effects in the central 
region do not affect its non-equilibrium statistics. 

Finally we would like to mention that it is however pos- 
sible to recover the Ng ansatz from our results in the limit 
of low-energy scales^, i.e. when (u> — /i Q ) ~ and then 
exp(w — fi a )/kT ~ 1. This implies that although approx- 
imate, the Ng ansatz might be good enough to describe 
low-energy excitations like, for example, the Kondo ef- 
fect in correlated electron systems, which gives a sharp 
feature in the spectral density around the Fermi level at 
equilibrium, or split Kondo peaks around the leads' Fermi 
levels at non equilibrium 70 . However, such an ansatz 
will most probably fail to describe systems for which 
the interaction (electron-phonon, electron-plasmon) is 
restricted on an energy scale defined by the phonon 
(plasmon) frequency, which is finite and not necessarily 
smalPEd 
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